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1. NSPT, Langevin equation, gauge fixing and all that 

Numerical Stochastic Perturbation Theory (for a review see Ref. [|l|]) is a powerful tool to study 
higher-loop contributions in Lattice Perturbation Theory (LPT). LPT is much more involved than 
perturbation theory in the continuum, and thus only few results beyond one-loop level are available. 
There have already been various applications of NSPT in the past: the average plaquette to very 
high orders in pure Yang-Mills theory to identify the gluon condensate [|2|], the residual mass for 
lattice HQEF renormalization factors for bilinear quark operators [Q], renormalization factors 
related to the QCD pressure etc. Relatively new is the application of NSPT to gluon and ghost 
propagators in Yang -Mills theory ^. Here we report on first steps towards an NSPT study of 
the ghost propagator in Landau gauge, in particular at two-loop level. 

It is known that the lattice Langevin equation with an additional running "time" t, beyond the 
four physical dimensions, leads to a distribution of the gauge link fields according to the measure 
exp (— ^oi?/]) in the limit t ^ oo. Discretizing the time t = nx and using the Euler scheme, the 
equation can be solved numerically by iteration: 

+ 1; T]) =exp(-F,,^[[/,T]]) U,_^{n■^) (1.1) 

with a force containing the gradient of So and a Gaussian random noise T] , 

Fv,^,[t/,r]] =i(TV,,^5G[t/] + ^/^T],,^). (1.2) 

Vjc,^ is the left Lie derivative acting on gauge group-valued variables while So is Wilson's one- 
plaquette gauge action. 

In NSPT one rescales £ = jSt and expands the link fields (and the force) in terms of the bare 
coupling constant g<x p^^/^: 

[/,,^ (? ; T] ) ^ 1 + £ ^ (? ; T] ) . (1.3) 

/>() 



Then the solution (1.1) transforms into a system of updates U U', one for each perturbative 
component U : 

[/(1)' = [/(1)-FW, [/(2)' = [/(2)_77(2) + ^(/7(1))2_^(1)^(1)^ ... (1.4) 

The random noise rj is fed in only through F^^\ higher orders become stochastic by propagation 
of noise through the fields of lower order. 

In terms of the (algebra- valued) gauge field variables A = log?/, 

A.,^,(?;T])-XiS"'^'^5(^;^)' 41 = rX'?, (1-5) 

/>0 

we are enforcing antihermiticity and tracelessness to all orders in g by requiring 

AWt = _^(')^ TrAW=0. (1.6) 

The Landau gauge is achieved by iterative gauge transformations using a perturbatively expanded 
version of the Fourier-accelerated gauge-fixing method |j8|] applied to each 50-th configuration 
in the Langevin process. Only these are evaluated in order to control the autocorrelations. Each 
Langevin update ( 1 .4 ) is completed by a stochastic gauge-fixing step and by subtracting zero modes 
of A^') as described in Ref. [Illl. 
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2. The ghost propagator in NSPT and in standard LPT 

The continuum ghost propagator G{q^) in momentum space is defined as G"^{q) = 8"^G{q^). 
On the lattice it is obtained as the color trace 

G{aq{k)) = ^±-^G'^"[aq{k)) = {TvM-\k))^ (2.1) 

as a function of the lattice momenta aq^{k) = Ink^a/L^ associated with plane waves |^) labelled 
by integers k^ = (— L^/2,L^/2] . In Landau gauge, the ghost propagator requires the computation 
of the inverse of the Faddeev-Popov (FP) operator 

M = -d-D{U), (2.2) 

with D{U) being the lattice covariant derivative and d the left lattice partial derivative. M^^{k) in 
( p] ) is the Fourier transform of the inverse FP operator. 
The perturbative expansion is based on the mapping 

{A^J^^}^{M^')}^{[M-'f}. (2.3) 

With an expansion of M in terms of M^'^ containing , a recursive inversion is possible in coor- 
dinate space: 



[M-'f^ = m(0) ' , [M-'Y'' = - [M^'Y' £m('-^') [M-'Y> . (2.4) 

j=0 

The momentum-space ghost propagator at n-loop order is obtained from even orders I = 2nof ' 
sandwiching its foregoing expansion between the plane- wave vectors: 

G^''\aq{k)) = {k\ [M-^f^^"^\k). (2.5) 

Odd / orders have to vanish numerically. We discuss the results in terms of two forms of the 
dressing function for one and two loops: 

j'^''\aq) = {aqfG'^"\aq{k)), J^"\q) = q^ G^"\aq{k)) . (2.6) 

Here we use the standard notation for hat-variables, e.g. 

In standard LPT, loop contributions are calculated in the infinite volume and a — > limit. In 
this limit the two dressing functions coincide. The renormalization of the dressing function is 
performed in the RI'-MOM scheme: 



with the renormalization condition 



J^\q,pi,aKV%2^^2 = \. (2.9) 
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Restricting ourselves to two-loop order, we have e.g. 

2 



' / / 1 

i=\ k=o \^ 



(aq)' 



(2.10) 



Only the leading coefficients zf]^' are entirely calculable in continuum perturbation theory (PT): 
Zj^f = —3Nc/2, Z22 = —35N^/S. The non-leading coefficients zfj^'\i>k>Q are only partly known 
from PT: z.-^( ~ + T^i^o) ' ^n) ^^'^^ to be calculated in LPT For example, entering Z2\ 

is z^Q = 13.8257, known from one-loop LPT [^], while z^q is unknown. 

From the relation [|^] Wri/ = tto + (-(22/3)A^flog(a/i) +73.9355) «(,+..., with the bare 
couphng tto = A'^c/ (8^^j8), we get for the two-loop dressing function: 

72-l°°P(a,^,^) = 1 + i (7i,ilogMH/l,o) +^ (/2,2log^M^+/2,llogK)^+/2,0)(2.11) 

with 

7ii = -0.0854897 , 7i,o = 0.525314, 72,2 =0.0215195, 72,i = -0.358423 (2.12) 
and the unknown finite two-loop finite constant /2,o or ^2^' , 

72,0 = L47572 + 0.00144365 zf^i . (2.13) 



3. Results 

The aim of this first investigation of the ghost propagator in NSPT was the confirmation of 
the known 7i_o, and a prediction of the unknown 72,o- We concentrate ourselves on an analysis of 

As an example of the measured ghost propagator we show the one- and two-loop results /(^^ 
and /(^) for the dressing function in Fig. |I] together with / ("=^/2) that is bound to vanish. 
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Figure 1: Measured ghost dressing function J{q) vs. q^ for all inequivalent lattice momentum 4-tuples 
(fci ,^2,^3,^4) - see (2.2) - near the diagonal ones for lattice sizes L = 6, . . . , 14 and for the time step e = 0.01. 
Left: The one-loop /3^') and two-loop (<=<: /3^^) contributions, right: the vanishing contribution jS^^/^. 
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3.1 The limits to be taken 



• The Umit £ ^ 0: We solved the Langevin equations for different step sizes e = 0.07, . . . , 0.01 
and obtained the Langevin result for each chosen momentum set of the propagator at fixed 
lattice size L and £ = by extrapolation as shown in Fig. ^. 




E E 

Figure 2: Linear plus quadratic correction extrapolation to e = of the one-loop (left) and two-loop (right) 
ghost dressing function for the momentum tuple (1,1, 1,1) on a lattice of size 12^. 



• The limits L^oo and a ^0: In order to make contact with standard LPT both Umits have to 
be performed. To extract the non-logarithmic constants in those limits we make the following 
ansatz for the dressing function taking into account hypercubic symmetry (one-loop example; 
here we use the standard notation for hypercubic invariants) 

J^'\q)=Ji.riogf"+Jlo-Li^), (3.1) 

^4 ^6 

JhO-dq) = J\fi- L + c\q^ + C2^+c^q^ +C4 {q^f + cs + ce {q^f H (3.2) 

q^ q^ 

The problem arising here is how to represent - on finite lattices - the logs that appear in the 
L^oo regime. Our proposal here is to replace the divergent lattice integrals, that give rise to 
the logarithms, by finite lattice sums and use these expressions in the fits at fixed L. 

3.2 Handling the lattice logs encountered 

We illustrate the procedure by the example of a typical one-loop divergent integral 

/7l/a aAi. 1 

In the limit aq ^0 [ pT] ] one gets 

A{aq) = -log{aq)^ + au ay = 2 + Fq- Je = 5.19201 . (3.4) 



On a lattice with finite L we calculate the corresponding lattice sums: 

1 



A(/^L) = 1 X ^ 



[Itisin^(f(/v-/v^))] 



(3.5) 
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with ak. 



I > {'V ) £ (~ f > f ] • This leads - for each L - to the replacement: 



7i,olog(a^)^ ^2 7i,o(A(/'?,L)-ai) . 



(3.6) 



This also results in a reshuffling of irrelevant terms. The result is a flattening of the data with the 
log-terms subtracted (see Fig. ^). This then allows to extract the V ^ o° limit fitting the remaining 
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Figure 3: Original and remaining "non-logarithmic" contributions to J using logarithms and lattice loga- 
rithms at one-loop and two-loop level as function of for a lattice lA"^. 



non-logarithmic data (at present no momentum cuts on the data are used) with the ansatz (3.2). In a 



similar spirit, a log-squared behavior in a two-loop contribution is modeled by using the following 



expression as a discretized version of [ ]llp 

r^l" d^k 1 



E(ac]) = (4ny -^-4i ^- TA(aA:)^-log2(a^)2-(ai + l)log(a^)2 + 28.0086 (3.7) 

J-nlai27ifp^y 2 

where A{ak) and are defined in ( p.3[ ) and (3.4) 



3.3 Results based on the outlined fitting procedure 

The results for /i,o;L and /2,o;l as function of l/L"* are shown in Fig. ^ A linear fit for 



L = 10, 12, 14 leads to the one-loop result 
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Figure 4: The V ^ °° limit of the constant J„,o-^l- 
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77^ = 0.5255(24) (3.8) 

in agreement with the expectations. A linear fit as in the one-loop case would lead to a preliminary 
two-loop value /j^q = 1.47(2). This results in the non-logarithmic contribution z^q to the two-loop 
ghost self-energy in the RF-MOM scheme in Landau gauge being compatible with zero. 

4. Summary 

• We have performed the first two-loop calculation of the lattice ghost propagator in Landau 
gauge. 

• The one-loop constant 7i o agrees with the known V ^ oo result. 

• The two-loop constant 72,o has been estimated for the first time. 

• A detailed analysis of all necessary limits has been performed. 

• A proposal about how to mimic the usual logarithmic terms on finite lattices is made. An 
alternative procedure outlined in Ref. is under development. 

• A detailed comparison for a finite volume and a set of lattice momenta with Monte Carlo 
data would be desirable in order to separate out the nonperturbative effects on the ghost 
propagator. 
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